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Abstract 

Diffraction radiation (DR) is produced when a charged particle passes through an 
aperture or near a discontinuity in the media in which it is traveling. DR is closely 
related to transition radiation (TR), which is produced when a charged particle tra- 
verses the boundary between media with different dielectric constants. In contrast to 
TR, which is now extensively used for beam diagnostic purposes, the potential of DR 
as a non-interceptive, multi-parameter beam diagnostic remains largely undeveloped. 
For diagnostic measurements it is useful to observe backward reflected DR from an cir- 
cular aperture or slit inclined with respect to the beam velocity. However, up to now, 
well founded equations for the spectral-angular intensities of backward DR from such 
apertures have not been available. We present a new derivation of the spectral angular 
intensity of backward DR produced from an inclined slit for two orientations of the 
slit axis, i.e. perpendicular and parallel to the plane of incidence. Our mathematical 
approach is generally applicable to any geometry and simpler than the Wiener Hopf 
method previously used to calculate DR from single edges. Our results for the slit are 
applied to the measurement of orthogonal beam size and divergence components. We 
discuss the problem of separating the simultaneous effects of these beam parameters on 
the angular distribution of DR and provide solutions to this difficulty. These include 
use of the horizontal and vertical polarization components of the radiation from a sin- 
gle slit and interferences from two inclined slits. Examples of DR diagnostics for a 500 
MeV beam are presented and the current limitations of the technique are discussed. 
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1 Introduction 



The new generation very high power radiation devices such as short wavelength free electron 
lasers and single pass spontaneous emission (SASE) devices for the generation of intense uv 
and x-ray beams requires the development of linear accelerators which produce optical quality 
charged particle beams. The high power, small size and low emittance of these beams present 
an enormous challenge for both the diagnostic measurement of beam parameters and the 
accurate positioning and control of these beams. Conventional screens or other interceptive 
probes are incompatible with the operation of such accelerators in many instances. Non 
interceptive devices such as wall monitor arrays, give limited beam information. Other 
types of monitors, such as synchrotron monitors, while non-invasive, cannot be used in 
linear beam line geometries. Hence the development of low cost, compact, nondestructive 
monitors capable of measuring multiple beam parameters would be very useful for many 
high power beam applications. 

We have investigated how the properties of diffraction radiation (DR) can be used to 
measure beam divergence, energy, position, transverse beam size and emittance. DR devices 
such as slits or circular apertures through which the beam passes offer minimal perturbation 
to the beam, respond rapidly to changes in beam parameters and are inherently compact, 
so that they can be implemented at many places in the beam line. 

Diffraction radiation (DR) is produced when a charged particle passes through an aper- 
ture or near an edge or interface between media with different dielectric properties. DR 
is closely related to transition radiation (TR), which is produced when a charged particle 
traverses the boundary between media with different dielectric constants. DR has been 
studied theoretically since the early 1960's, [1-5]. However, in contrast to TR, which has 
been well investigated experimentally and is now used extensively for beam diagnosis, we are 
aware of only two experimental studies of DR in the literature [6], [7]. In only one of these 
investigations [6] was DR actually used as a beam parameter diagnostic, i.e. to measure 
the beam bunch length. Furthermore, in this study forward and 180° backward directed 
DR and TR were observed simultaneously. This situation complicated the comparison of 
measured to theoretically predicted DR properties. Therefore, to avoid such difficulties and 
for experimental convenience, it is advantageous to observe the radiation from an inclined 
aperture, e.g. a slit inclined at 45°, which produces DR both in the forward direction and 
in the direction of specular reflection of the virtual photons associated with the particle, i.e. 
at 90 degrees with respect to the beam velocity. 

The detailed nature of the far field angular distribution (AD) of forward directed DR 
has been derived and discussed by Ter-Mikaelian [4] for an electron normally incident on a 
circular aperture and a slit in an infinite, thin plane. The fields and AD of DR from an 
electron incident on a single edge (half plane) inclined at an arbitrary angle with respect to 
the velocity vector of the particle have also been derived and discussed in Refs. [2], [4] and [8]. 
Although some authors have surmised that reflected DR from an inclined circular aperture 
or slit has properties similar to reflected TR [9-10], a derivation, from first principles, of the 
spectral-angular distribution of reflected DR from such apertures has not been previously 
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available in the literature. We present here solutions for the fields and the horizontally and 
vertically polarized intensities of DR from an inclined slit for two orientations of the slit 
edge, i.e. parallel and perpendicular with respect to the plane of incidence. 

Previous theoretical investigations of forward directed DR generated from electron beams 
passing through circular apertures and slits [9-12] have examined the potential use of the far 
field angular distribution (AD) to diagnose beam properties such as position with respect to 
the center of the aperture, beam size, trajectory angle, divergence and energy. The potential 
advantage and difficulty in using the AD of DR for beam diagnostics is that it depends on 
the beam position (size) as well as the beam divergence and energy. In contrast, the AD 
of TR depends only on the beam divergence and energy. The challenge for the application 
of the AD of DR as a diagnostic then is to devise a means to separate out the competing 
effects of divergence and beam size on this observable. 

We have developed a number of ways to address and solve this problem. We demon- 
strate here how the measurements of the horizontal and vertical polarization components 
of backward reflected DR from a single slit can be used to separate out the effects of beam 
size perpendicular to the slit edge and the component of divergence parallel to the slit edge. 
Furthermore, we show how the measured AD's for two different orthogonal orientations of 
the slit axis, can be used to determine the components of beam size and divergence in two 
orthogonal directions, i.e. the directions perpendicular and parallel to the plane of incidence. 
In addition we suggest how the interferences from two slits can be used to make a sensitive 
measurement of divergence alone. 

2 Diffraction Radiation Properties 

When a charged particle or ensemble of particles passes through an aperture or near an edge 
at a distance d, appreciable DR is generated when the condition: d < 7 A is fulfilled, where 
A = A/27T, and A is the observed wavelength. Table I. presents typical wavelengths of DR 
observed from a slit whose width a = 2d = 2mm for a variety of lepton and hadron beams. 
Table I. shows, for example, that DR is produced in the optical part of the spectrum (i.e. 
IR to visible) for electrons with energies in the range of 0.5 - 5 GeV. 

Table I. 

Typical wavelengths of diffraction radiation generated from a 2-mm wide slit 

for various accelerators 



E(GeV) 


7 


A(/xm) 


Accelerator 


0.5 


10 3 


6.30 


APS 


5.0 


10 4 


0.630 


CEBA 


50.0 


10 5 


0.063 


SLAC 


250.0 


250 


25.00 


RHIC 


10 3 


10 3 


6.30 


FNAL 
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When the observed DR wavelength is much smaller than the transverse or longitudinal 
size of a charged particle bunch, incoherent DR is produced independently from each particle 
and the total intensity / ~ N, the number of particles in the bunch. For wavelengths 
comparable to the bunch dimensions coherent DR is observed with an intensity / ~ N 2 . The 
DR diagnostic methods we will present can utilize either incoherent or coherent DR. However, 
for simplicity and to illustrate and emphasize the basic methodology we will concern ourselves 
here only with the analysis of incoherent DR. Furthermore, since we have discussed DR 
diagnostics using circular apertures elsewhere [10], we will focus our attention in this paper 
to the generation of DR from slits. An extended analysis of coherent DR from a slit will be 
presented in a future work. 

3 Theory of Backward Diffraction Radiation from an 
Inclined Slit 

3.1 Parallel orientation of slit with respect to the plane of incidence 

A relativistic charged particle passing through a slit in a conducting screen emits diffraction 
radiation in the forward and backward directions, the latter in the direction of specular 
reflection of the virtual photons associated with the relativistic particle. The peak intensity 
of diffraction radiation, like transition radiation, occurs at an angle 9 ~ 7 -1 , where 9 is 
measured from the direction of the particle's velocity vector v for forward DR, and from the 
direction of specular reflection for backward DR. 

Figure 1A. shows a side view, and Figure IB. a top view, of the cone of backward 
diffraction radiation emitted by a charge passing through a slit which is inclined at the angle 
\l/ with respect to the velocity vector and oriented such that its edge is horizontal or parallel 
to the plane of incidence (n, v plane), where n is the unit vector normal to the plane of the 
slit, i.e. the (x 1 , y') plane. Additionally, the particle passes through the slit with an offset 5 in 
the vertical {%/') direction measured from the center of the slit which is defined by the x' axis. 
For simplicity, we have not drawn the forward directed cone of DR which is also generated in 
the direction v. Figures 1A. and IB. show the coordinate system (x,y,z) used to describe 
the radiation fields associated with the wave vector k observed in the backward direction 
(i.e. directed along the z axis). The slit is assumed to be an infinitely thin conducting screen 
with infinite extent in the x' and y' directions beyond the slit aperture. 

Our calculation of the backward reflected DR fields employs three different concepts (a) 
the Huygens-Fresnel diffraction approach, following the treatment given by Ter-Mikaelian [4] 
for the case of forward DR, (b) the concept of backward reflected pseudo photons, used by 
Wartski [13] to calculate reflected transition radiation and (c) a form of Babinet's principle. 
Our approach is quite different and mathematically much less complex than the exact Wiener 
Hopf method used by previous researchers to calculate backward DR from a single edge [2,4]. 
In addition, our method is generally applicable to other geometries where the Wiener Hopf 
technique is very difficult to implement. 
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Consider a thin reflecting metal screen ^with infinite extent. A charged particle crossing 
this surface produces the well known forward and backward transition radiation. Compare 
this to the case of Fig. 1., where the slit has an area Si and the conducting screen has an 
area S 2 . Babinet's principle states that the transition radiation (TR) field from the infinite 
conducting screen, S^, can be obtained from a Huygen's Fresnel integral over the infinite 
area S^ — Si + S 2 , so that the TR field can be expressed asE 00 =Ei+E 2 where E x 
and E 2 are those parts of the field which are obtained from a Huygen's Fresnel integral 
over the areas Si and S 2 , respectively [8]. If material of area Si is removed, then the field 
E 2 = E oo — E 1, where E 2 is the desired backward DR field. Note that if the area Si — > 0, 
E 2 — > E oo , and the DR field approaches that of TR. 

The incident field in the Huygens-Fresnel integral is the field associated with the rela- 
tivistic particle, Ei . Using Ter-Mikaelian's approach to calculate the diffraction radiation 
field and Wartski's method of calculating backward reflected radiation, we can express the 
x and y components of the radiation field as 

E X:y (k x ,k y ) = ni,±(w,*)^2 J j E iXiiy (x,y)e- ik -~?dxdy (1) 

where E ix>iy are the components of the field of the incoming charged particle and k x and k y are 
the components of the wave vector k in the plane normal to the direction z shown in Figure 
1A. and IB. As is the case for TR, there are generally three components to the radiation: 
direct, reflected and transmitted. However, for an opaque, highly conducting screen, the 
dominant term for backward TR or DR production is the reflected component. The fields E x 
and E y for backward TR and DR are thus proportional to the Fresnel reflection coefficients 
r\\ and r±, respectively. In Eq. (1), the integral over the area S 2 is done in terms of the 
coordinates (x, y) of Figures 1 and 2. The far field reflected spectral intensity of DR is found 
by integrating the flux of pseudo photons reflected off the surface S 2 , which pass through 
the plane that is perpendicular to the direction z, i.e. the direction of specular reflection of 
incident pseudo photons having a wave vector k, which is parallel or closely parallel to the 
velocity vector v. The phase term in Eq. (1), which is a function of ~~p = (x',y'), a vector 
lying in the plane of the screen, can be related to the coordinates (x, y) by the transformations 
x = x sin(\I/), y = y , and the differential elements of area by dxdy = dx dy sin(\E'). The 
detailed solution of the integral in Eq. (1) is presented in the Appendix. 
The results are: 



4tt 2 v f 



_ e ~ a i(J~ ik y) _| \ g-"2(/+ifc w ) 



(2) 



{f -ik y ) (f + iky) 

which is identical to Eq. (31.18) of Ref. [4], except for the reflection coefficient r^o;,^), 
and the scalar quantity k x = k x — \k/{2^ 2 )] cot and 



E (k k ) — e l±tl3 



_ e - a i(f- ik v) e - a '2(f+ik y ) 



(/ - iky) (/ + iky) 



(3) 
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which is identical to Eq. (31.19) of Ref. [4], except for the factor r±(uj, ty). Here a ii2 = 

a/2±5, f = [k x + uj 2 / (v 2 ^ 2 )} 1 ^ 2 , v = |v|, k x = ksm.9 x , k y = ksm9 y and 9 XjV are the projected 
angles of the vector k into the x, z and y, z planes, respectively. 

The horizontally polarized intensity (i.e. parallel to the slit edge) and vertically polarized 
intensity (i.e. perpendicular to the slit edge), which are observed in the plane perpendicular 
to the direction of specular reflection, are defined in terms of the x and y components of the 
fields: 

d 2 N horiz . 2n 2 k 2 c 2 

I As | (4) 



dwdfl fiw 

I |2 



ak 2 k 2 x e- af 
A^uoP{P + kl) 



[cosh(2/£) + sm(ak y + &(k x , k y ))} 



and 



d 2 N vert . 2n 2 k 2 c 
dudfl hjj 



respectively, where 



Ey\ Z (5) 



|2 ak 2 e~ a f 



' r± ' Att 2 uj (p + k 2 ) [ cosh ( 2 ^) ~ sin ( ak y + H k x, k y ))} 



k y ) = sm-^f - kl)/(f 2 + ^)] = co^[-2fk y /{f + k% (6) 



and a ~ is the fine structure constant. 



The equations for vertically polarized forward diffraction radiation intensity, first ob- 
tained by Ter-Mikaelian [4] and later rewritten and presented in Ref. [9], are both equiva- 
lent to Eq. (5 ) without the factor |rj_| 2 . Note, however, that the equations for phase term, 
&(k x ,k y ) presented in Ref. [9] are incorrect. The correct relationships are given in Eq. (6). 



3.1.1 Limiting Forms for Backward DR at 45° 

Under the usual small angle approximation for observation of TR and DR from relativistic 
particles, i.e. v m c,9 ~ 7~ x 1, the quantities k x k x k9 X) k y ^ k9 y , and k z ^ k. We 
then can consider several limiting forms for backward DR. 

(1) Single Edge Limit: When either a\ or a 2 — > oo , the sinusoidal term in Eqs. (4) 
and (5), which represents the interference between the intensities from each edge of the slit, 
disappears and expressions for the horizontal and vertical intensities of DR from a single 
edge are obtained. The addition of the horizontal and vertical intensities taken in this limit 
produces an expression for the total DR intensity which is identical to that obtained using 
the exact Wiener Hopf approach (see Ref. [8], Eq. (14)). The correspondence of the results 
of our calculations to those calculated using exact theory is a strong confirmation of the 
validity of our method for calculating backward reflected DR. 
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(2) Small Displacement Limit: When 5, the particle displacement from the center of the 
slit, is small, i.e. 5 <C 7 A ~ a, the hyperbolic cosine in Eqs. (4) and (5) can be expanded 
in a Taylor series. Retaining terms only up to second order, and using the small angle 
approximations given above, Eqs. (4) and (5) become: 

d 2 N horiz . = . |2 _a_ a X 2 e -m+W 

dcudn " |r|11 4tt 2 w 7 (i + x 2 )(i + x 2 + r 2 ) u 

•[1 + 2(^) 2 (1 + X 2 ) + sm(RY + $'(X, Y))] 
7A 

and 



d 2 N ♦ o a e -R(i+x*y/* 

u ly vert. I 1 2 Lt 2 /oN 

dwdn |r±l 4^J 7 (1 + X 2 + F 2 ) 1 J 

•[1 + 2( X) 2 (l + X 2 ) - sm(RY + $'(X, Y))} 
7A 

where we have introduced the new variables X = ^Q x and Y = 7^, the x and y projected 
angles scaled in units of 7 _1 , R = a/^\ and the reduced phase term 

&'(X,Y)) = sin- 1 [(l + X 2 -F 2 )/(l+X 2 + F 2 )] = cos- 1 [-2(l+X 2 ) 1 / 2 F/(l + X 2 + F 2 )]. (9) 

(3) TR Limit: Additionally, when R « 1 and 5 « 7A , Eqs. (7) and (8) each reduce to 
one half the intensity of TR, which is the expected correct limit. In this regime the particle 
radiates as if the slit in the screen were absent. 



3.2 Perpendicular orientation of the slit with respect to the plane of incidence 

Figures 2A. and 2B. illustrate the case in which the slit is oriented such that the edge is 
perpendicular to the plane of incidence, i.e. the plane containing v and n the normal to the 
plane of the screen. In this case, the offset e of the particle velocity vector from the center 
of the slit is in the direction Y' and the radiation fields E x and E y can be shown to be the 
same as those given above in Eqs. (2) and (3) but with a\ — > aisin^, 02 - * 02 sin ^ and 
5 — > esin 1 !/. Also, E x oc r±(uj,^) and E y oc ri|(a;,\I f ) since, for this orientation of the slit, 

E x is perpendicular and E y is parallel, respectively, to the plane of incidence. Then E x and 

E y take the same forms as Eqs. (2) and (3), with the coefficients r± and r\\ interchanged. 

Similarly the vertical and horizontal intensity components take on the same forms as Eqs. 

2 1 1 2 

(4) and (5), with the coefficients |r_i_| and ry interchanged. 

In Ref. [8] it has been shown that an exact calculation of reflected DR produced by a 
particle traveling with a velocity with a parallel component along the direction of the edge 
of a semi-infinite screen [2] reduces, in the limit 7 >> 1, to the same result as that of a 
particle whose velocity vector is purely perpendicular to the edge (see Eqs. (7) and (14) of 
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Ref. [8]). The former situation corresponds to the geometry of Figure 1. as discussed above in 
Section 3.1; the latter situation corresponds to the geometry of Figure 2., which is discussed 
in this Section. As in the parallel case, the single edge results for the geometry of Fig. 2. 
are completely reproduced using our approach in the limit ai or a2 — > oo. Furthermore, in 
the perpendicular incidence case the small angle, relativistic limit produces the same forms 
for the horizontal and vertical intensities as the parallel incidence case, i.e. Eqs. (7) and 

2 I 1 2 

(8), but with the coefficients |rj_| and rJ interchanged. In addition, when R « 1, the 
horizontal and vertical intensities each go to the same TR limit as in the parallel incidence 
case. These results provide further evidence of the soundness of our theoretical approach. 

For the sake of simplicity and to illustrate the main results of our analysis of DR from 
slits, we will present numerical examples only for the case of the slit oriented with its edge 
parallel to the plane of incidence . However, numerical results can easily be obtained for the 
perpendicular orientation by simple substitution of variables as listed above. 

4 Discussion and Computational Results 
4.1 General properties of DR 

Eqs. (5) and (6) indicate that the highest yield of backward DR observed in the IR to visible 
part of the spectrum will be generated when the surface of the slit is a highly conducting 
mirrored surface. In the case |r_i_| — \r\\\ ~ 1. and for 7 >> 1, Eqs. (7) and (8) show 
that the polarization components are the same as forward DR. These equations also reveal 
some interesting characteristics of DR from a slit, which have not been fully discussed in 
previous studies, and which can be easily visualized with the help of three dimensional plots 
of intensity observed in the X, Y plane. 

Figs. 3 A. and 3B. show the horizontal and vertical components, respectively, of the 
intensity of DR for a value of the ratio R = a/ (7A) = 0.5. Similarly, Figs. 4A. and 4B. show 
the component intensities for R = 2.0 . The intensity axis (Z) in these Figures is scaled in 
proportion to the corresponding intensity component of transition radiation. For example, 
Fig. 3A indicates that the peak intensities of horizontal and vertical DR, respectively, are 
about 25% and 80% that of TR , when R = 0.5. Figures 3. and 4. also show that diffraction 
fringes are present in both polarization components, but are only observed in the Y direction 
(perpendicular to the slit edge) for both the orientations shown in Figs. 1 and 2. This can be 
deduced directly by an examination of the term sin(RY + &'(X, Y)) in Eqs. (7) and (8), i.e. 
the frequency of oscillation is determined solely by the value of k y component of the wave 
vector represented through the variable R. Secondly, the inteferences in the horizontal and 
vertical intensities are 180° out of phase. This is evident by the presence of the difference 
in the sign of the sinusoidal terms in each component. Thus, in general, the maxima of the 
vertical component along the Y axis occur at the minima of the horizontal component. 

One further observes that the equation for horizontal intensity, Eq.(7) contains the factor 
X 2 /(l + X 2 ) , which is not present in Eq.(8) for the vertical intensity. This term, which 
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is similar to the form of TR, forces the horizontal component to have a null at X = for 
every value of Y . A comparison of the ratio of vertical to horizontal intensities reveals that 
the vertical component is larger than the horizontal component, and that ratio becomes 
smaller as R decreases approaching unity as R — > , i.e. the transition radiation limit. 
Also, as R becomes large the DR intensity drops off exponentially. Thus the value of R 
should be kept to a minimum to maximize the observed intensity of DR. We have chosen 
the values R = 0.5 and 2.0 to show that measurable horizontal and vertical intensity levels 
can be achieved for reasonable experimental conditions and to demonstrate that the ratio 
of vertical to horizontal intensity depends on the value of R. For example, a value of R = 2 
applies to the case when the beam energy E b = 500 MeV, the slit width a = 2mm and the 
wavelength A ~ 3/im. By change of the parameters a and A, or by use of different values 
of R , other workable combinations are possible for a given beam energy. Note that as the 
beam energy increases, a smaller value of R can be achieved for the same wavelength and 
slit size. Since the intensity of DR, like TR, is proportional to 7 2 (see Eqs. (7) and (8)), DR 
becomes progressively more intense as hence easier to detect as the beam energy increases. 

In general, the sum of horizontal and vertical polarization intensity components will 
be observed by a detector in the X, Y plane. Thus, if the angular distribution of DR is 
imaged, e.g. by using a camera with a lens focussed at infinity, the sum of horizontal and 
vertical distributions will be superimposed. However, either distribution can be imaged 
separately by placing a rotatable polarizer in front of the lens. We will show that each 
of these intensity components provides information about the beam size and divergence in 
each of two orthogonal directions and, furthermore, that all of these beam parameters can 
be measured by proper analysis of the horizontally and vertically polarized intensities. 

4.2 Effect of Beam Parameters on the Angular Distribution of DR 
4.2.1 Beam Size 

The above discussion shows that both the horizontal and vertical DR spectral-angular dis- 
tributions are functions of the beam position 5 relative to the center of the slit. This effect 
has been previously considered [9,11] as a possible diagnostic of the beam position and size. 
In Ref. [9] it was shown that the effects of beam size and offset on the angular distribution 
of DR are the same. Also, in Ref. [11] it was shown that the DR intensity is a minimum for 
a beam centered in the slit, i.e. zero offset. Since the beam can be centered by monitoring 
and steering the beam to minimize the total intensity, the offset can be nullified or otherwise 
determined, for example, by use of a standard wall current monitor. The remaining effects on 
the angular distribution of the DR will then be the beam size in the direction perpendicular 
to the slit edge and the two orthogonal components of the beam divergence. 

To determine the effect of beam size on the intensity components, we assume that the 
beam centroid has been centered in the slit by means of one of the previous methods described 
above and that any particle within the finite beam spatial distribution has the same effect 
as the offset of a single particle. We further assume that the beam has a separable spatial 
distribution given by S = Si (5) ■ S^e), where 5 and e are in the directions y' and x', 
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respectively. For the parallel orientation of the slit, 5 and e are perpendicular and parallel, 
respectively, to the slit edge. Integration of the distribution S over 8 and e produces the 
average beam sizes (5) and (e). Because of the simple dependence of the above expressions 
for horizontal and vertical intensity on the variable 6, an integration of S over the intensity 
components is mathematically equivalent to a simple replacement the variable 5 with its 
average value (5) . To simplify the notation in the remainder of the text we will identify and 
refer to S and e as the othogonal rms beam size components, which are perpendicular and 
parallel, respectively, to the plane of incidence. 

4.2.2 Beam Divergence 

Ref. [9] has examined the effect of beam size perpendicular to the edge of a slit on the vertical 
intensity component of forward DR in the absence of divergence. A discussion of the effect 
of beam size on the horizontal component was dismissed because of its much lower intensity 
in comparison to the vertical intensity - an effect which was due to the rather large value of 
R = 2ir used in the analysis. However, for any realistic beam, the divergence is not negligible 
and affects both the horizontal and vertical intensities of DR. Furthermore, in general, the 
horizontal component of DR is not negligible in comparison to the vertical component and, 
as we will show, this component provides important information on the component of beam 
divergence parallel to the slit edge. It is therefore important to analyze both the effect of 
beam size and divergence on the vertical and horizontal intensity components to fully assess 
the diagnostic potential of DR. 

The effect of divergence must be taken into account by performing a two dimensional 
convolution of a distribution of particle trajectory angles projected in the X, Z and Y, Z 
planes (e.g. separable Gaussian distributions in 9 X and 9 y ). However, we have shown by 
numerical calculation that the effect of the rms divergence 8' on the horizontal component 
of the DR intensity, and the effect of the rms divergence e' on the vertical component is 
insignificant for S',e' < 0.2, a value which is large by the standards of most high quality 
accelerators. Then, to a good approximation, a one dimensional convolution of a line scan 
of the horizontal or vertical intensity calculated either in the plane defined by Y = const. 
or X = const, can in principal be used to predict the effects of the beam divergences e' 
and 5', respectively. Inversely, these divergences can be determined by fitting the convolved 
intensities to measured data. 

Horizontal Intensity Component The effect of the divergence e', the component par- 
allel to the incidence plane, on the horizontal intensity component can be estimated by 
neglecting the beam size term. Consider the term X 2 in the numerator of the horizontal 
component Eq.(7). Because of the presence of this term, a convolution of the horizontal 
intensity component with a distribution of particle angles G(X, e') will have the maximum 
effect on this component. Furthermore, the effect of e' will be maximized near the cen- 
ter of the pattern where the horizontal intensity goes to zero in the absence of divergence. 
One can estimate this effect by setting the variable Y in Eq. (7) equal to a constant value 
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Y — 5' — const., and noting the effect on a line scan of the one dimensional convolution of 
the horizontal intensity with e.g. the Gaussian distribution: G(X,e') = , 1 , 2 exp( 



-x 2 

2e> 2 >■ 

Figure 5. shows two such normalized scans of a ID convolution of horizontal intensity 
(Eq. (7)) with G for two different values of the variable, Y — 5' — 0.2 and Y — 0. In 
this calculation, R = 2.0 and e' = 0.1. As Fig. 5. shows, there is no discernible difference 
between the two patterns. Then only a small error will be made by setting the variable Y 
equal to a constant and performing a one dimensional convolution over X to infer the effect 
of the divergence e'. Inversely, it is possible to measure the divergence, e' from line scans of 
horizontally polarized DR for a number of Y values by comparing the measured intensity 
scans to theoretically convolved line scans in the variable X , which are parameterized by 
the variable e' . 

Now consider the effect of the beam size 5 measured perpendicular the plane of incidence, 
which is also perpendicular to the slit edge for the orientation shown in Fig. 1., on the 
horizontal intensity. Note that the effect of beam size e or any displacement parallel to the 
slit edge has no effect on this component. To compare the effect of the beam size S with effect 
of the divergence e', consider a line scan of the horizontal intensity component observed over 
a finite bandwidth (Au;) taken in the plane, Y — 0. Eq. (7) then reduces to: 

dN ^z. = I 12 «f AW X 2 R{1+X 2^2 _5_ 2 2 

dn 1 111 2vr 2 u (1 + X 2 ) 2 [ +l 7 A J 1 )V 1 } 

A convolution of this expression with the Gaussian distribution, G(X,e'), evaluated at 
X = 0, to first order in e' 2 , gives 

^^®G(X,.')oc.' 2 [l + (-i) 2 ]. (11) 

Thus, while it is still present, the effect of the beam size S on the horizontal component 
will be much smaller than that of the divergence e'. Figure 6. shows how a change in the 
divergence e' affects the horizontal intensity for a fixed value of the beam size 5 = 300/i. As 
predicted from Eq.(ll), the effect of e' is maximum near the origin. Conversely, the shape 
of pattern shape in the vicinity of the origin can be used to separate out and measure the 
divergence e'. 

Figure 7. shows the effect of the beam size, 5 on a line scan of the horizontal intensity 
(Eq. (12)) taken in the X, Z plane, i.e. Y — 0, when e 1 — and R = 2. Figure 7. shows 
that a change in 5 has no effect on the horizontal component in the region near X = 0, but 
does affect the peak value of the intensity observed at \X\ ks 1 and the fall off of intensity 
for \X\ > 1, for 8 > 50//. Numerical calculations indicate similar variation of the intensities 
with beam size and divergence for R = 0.5. Therefore, for simplicity, we will only present 
below numerical results for a single value of the parameter R = 2. 

Vertical Intensity Component A procedure similar to the one described above can be 
used to show that the effect of the divergence e' has a minimal effect on the vertical intensity 
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distribution. Figure 8. shows two scans of the vertical intensity observed in Figure 4B in 
the planes X — e' — 0.2 and X = 0. As similarly shown in Fig. 5., the difference in the two 
scans is very small. Thus, it should be possible to use a one dimensional convolution to infer 
the effect of the divergence 5', on line scans of vertically polarized DR (measured in the Y 
direction) for a number of X values. 

By comparing the measured intensity scans to a set of one dimensional convolutions in 
Y parameterized by the divergence component 5', one would hope to be able to measure 
this quantity. However, for the vertical intensity component the effect of the beam size 5 is 
comparable to that of the divergence 5'. To see this consider an observation of the vertical 
intensity component (Eq. (8)) in the Y, Z plane i.e. X = 0. For a finite bandwidth (Alu) 
measurement, Eq. (8) produces the vertical angular distribution 



dNvprt | |2 a Auo 2 e R 

\r±\ -— 2 —l 2 ^ : tpt, (12) 



'Vert. I , 2 _«_£^ 2 

dQ |r±l 4vr 2 to 7 (1 + F 2 ) 



,5,2 2F sm( RY) — (1 — Y 2 ) cos( RY) . 

•i 1 + 2 ^> 2 + — + n - '• 

To compute the effect of the divergence 5' one must take the convolution of Eq. (12) with 
a distribution over the variable Y, e.g. the Gaussian distribution. G(Y, 5') = exp( : =^ 3 -). 
We have performed this computation numerically to see how the effects of beam size and 
divergence compare. Figure 9. shows a vertical intensity scan taken at X = 0, for a fixed 
value of the beam size 5 =200/i and R = 2.0 for several values of divergence 5'. Figure 10. 
shows the effect of a change in the beam size, 5 on similar scans for the divergence value 
5' = 0. A comparison of Figs. 9. and 10. clearly indicates that the effect of divergences 
5' < 0.2 will compete significantly with the effect of beam sizes 5 < 300;U on the angular 
distribution of the vertical intensity component. Thus, in general, neither effect can be 
neglected and the two effects are not separable for this component. 

4.3 Strategies for separating beam size and divergence effects 

4.3.1 Use of slit oriented perpendicular to the plane of incidence 

In order to help separate out the competing effects of beam size and divergence, we have 
considered two alternative strategies. One method is to rotate the slit (or insert another 
into the beam line) so that the slit edge is perpendicular to the plane of incidence (see Fig. 
2A.). In this configuration the horizontal intensity component of the DR, convolved with 
the Gaussian function, G(X,S') = ex P(^r), can be used to measure the divergence 

5', which is now the component of divergence measured parallel to the slit edge. For this 
orientation the horizontal component is highly sensitive to 5' and only weakly depend on the 
beam size e , which is now the component of the beam size perpendicular to the slit edge. 
Evaluating the convolved horizontal intensity at X = 0, we obtain 

^|p®G(X,5')oc5' 2 [l + (^) 2 ]. (13) 
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Eq.(13) is the direct analog of Eq.(ll) above which, as shown in the previous section, can be 
used to determine e'. For the perpendicular orientation of the slit edge, however, the center 
of the horizontal AD pattern is sensitive to, and can be used to measure, the divergence 
component 5'. 

The vertical DR intensity component for this orientation of the slit is: 



dN V ert. I |2 OL Au 2 t R 



dn 1 111 4tt 2 lu 1 (i + r 2 ) 

. e 2Y sm( RY) — (1 — Y 2 ) cos( RY) . 



(14) 



which is the direct analog of Eq. (12). A one dimensional convolution of Eq. (14) with 
the Gaussian distribution G(Y,e') = ^-^2 ex P(^r)^ produces a pattern that is sensitive to 
the divergence e' as well as the beam size e. Thus for the perpendicular slit orientation the 
divergence e' and beam size e have comparable effects and it is not possible to distinguish 
between them using the vertical intensity component alone. 

4.3.2 DR Interferometry 

The second strategy we have devised to separate out and measure divergence is to use the 
interference of DR produced from two slits inclined at 45° with respect to the beam velocity 
in a configuration which is the direct DR analogy to a Wartski OTR interferometer [13]. 
In such a system forward DR from the first slit reflects from the second slit surface and 
interferes with backward DR generated from the second slit. The interferences produced by 
the two intensities will be superimposed on the single slit DR intensity. Figure 11. shows 
one possible configuration for a DR interferometer composed of two slits oriented with their 
edges parallel to the plane of incidence. The equation for interference DR, in analogy to 
TR, is obtained by multiplying Eqs. (7) and (8) by an additional interference term which is 
due to the difference in phase between forward and backward DR. The expressions for the 
horizontal and vertical intensity components are of the form: 

dN^ dN^ T 

Uiv Horiz,Vert A uly Horiz,Vert ■ It ^ \ /i r\ 

" 4 3f; sm (777-) ( 15 ) 



dtt dtt y 2L 



v 



where the superscript / refers to the two slit (interferometer), and S to the single slit angular 
intensity distributions, respectively, Ly = X/tt^^ 2 + 9 2 ) is the coherence length in vacuum 
for TR and DR, which respresents the distance over which the particle's field and the TR 
or DR photon differ in phase by one radian and L is the path length between the slits. 

For a fixed wavelength and energy the single slit term varies slowly with the angle 6. Then 
a convolution of Eq.(15) with a distribution of beam angles will chiefly effect the interference 



term sin 2 (^y^). The interference fringe visibility is then a function of the angular divergence 



and is independent of beam size effects. Thus it should be possible to use DR interferences 
to obtain a measurement of the beam divergence alone. Also as is the case with interference 
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OTR, measurements of the two polarized components of interference DR can be used to 
measure orthogonal components of the beam divergence [14]. The interference fringes, which 
are superimposed on the single slit angular distribution, provide increased sensitivity to 
angular differences in the particle trajectory angles and therefore smaller values of divergence 
(5',e' < 0.05) can be measured with an interferometer than with a single slit alone. 

Figures 12. and 13. show the vertical and horizontal components of intensity, respectively, 
from a DR interferometer produced by an electron beam with energy E = 500 MeV, R = 2, 
A=3.2/im, <5=200/x and 5' =0.05. The separation distance L = 6Ly = 3 meters for this beam 
energy. In any DR interferometer the angular distribution of the forward DR from the first 
foil will be partially cut off by the presence of the second slit in such a configuration. If the 
path length between the slits is large, the total angular field will only be cut off in the center 
of the pattern by a small amount. Thus only a small fraction of the interference patterns 
(I/67 for the parameter range specified above) will be cut out. This effect is not shown in 
Figs. 12. and 13. 

In contrast to the very weak effect of a small value of divergence, i.e. 8' = 0.05, on single 
slit DR (see Figs. 7. and 9.), the effect of this divergence on the DR interference fringe 
visibility is clearly visible. Since the fringes are well modulated out to Y ~ 4, the cutoff of 
the field of view due to the second aperture will have a negligible effect on a measurement 
of the divergence, which uses the interference pattern. With the help of a second camera or 
imager, a simultaneous measurement of backward DR from the first slit can be made and 
thereby provide addition single slit data for extracting the beam size. 

At the present time the use of the angular distribution of DR for diagnostics is practically 
limited to beams with moderate to high energies by the vacuum coherence length Ly ~ 7 2 A. 
The DR produced by upstream sources such as beam line discontinuities or other apertures 
will destructively interfere with backward DR generated from a diagnostic aperture when 
the distance between the upstream source and the aperture is much less than Ly. Also, the 
successful application of interference DR to measure beam divergences 5',e' < 0.05 requires 
that the inter-aperture distance L > 5Ly. Therefore, the diagnostic techniques described in 
this paper are limited to observation wavelengths and beam energies which do not give rise 
to an impracticably long coherence length for the accelerator facility being used. The means 
to overcome this limit must be developed before the methods described here can be applied 
to beams with very high energy. 

5 Summary 

We have developed a new method to calculate backward reflected diffraction radiation from 
any type of aperture. Using this method we have derived the equations for backward hori- 
zontal and vertically polarized spectral-angular intensities of DR from a slit inclined at an 
arbitrary angle with respect to the particle velocity. We have obtained results for two orthog- 
onal orientations of the slit edges, i.e. parallel and perpendicular to the plane of incidence. 
The results of our calculations reduce in the single edge limit precisely to those previously 
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derived using the exact Wiener Hopf method. However, in contrast to the Wiener Hopf 
approach, which has been successfully used to calculate DR from one type of aperture only, 
i.e. the single edge, our approach is readily applicable to any type of aperture and is much 
simpler to employ. 

From an analysis of backward DR for the two orientations of the slit mentioned above, 
we have developed a number of strategies to determine four unknown beam parameters of 
interest: the two orthogonal beam divergences, 5', e', which are perpendicular and parallel 
to the incidence plane, respectively, and the two corresponding orthogonal components of 
the beam size: 8, e. These strategies require the measurement of the horizontally and ver- 
tically polarized intensities from a single slit or a double slit DR interferometer. We have 
demonstrated that one dimensional convolved scans of the horizontal intensities in either 
the X = const, or Y = const, plane can be used to separately measure 5' and e', respec- 
tively. With the divergences known, the corresponding beam size components S and e can 
be inferred. Since each of these scans can be produced at multiple angles in the plane of 
observation, a large amount of data is available which can be used to reduce the error in the 
measurement of all the beam parameters of interest. The stategies we have developed are 
useful for moderate to high energy lepton or hadron beams. 



6 Appendix 

To proceed to evaluate the integral presented in Eq. (1) we express E ix (x,y) by its Fourier 
transform, so that we have 

K(k x , k y ) = 1 1 ( ^ + ^ + _ 2) e^-^d^dxdy. (16) 

Here, a = l/(P r yX), and the phase term appearing in Eq. (16) above is 

(A? _ t) . y = ( k ' x - k x )x'dnV + (k' y - ky)y' (17) 

where _^ 

k x = k x + (k z - k' ■ — ) cot \P. (18) 

v 

Note that k' • if = u , where u is the frequency of the Fourier component of the field. 
This phase term takes into account the variation in phase with position ~~p of the k! Fourier 
component of the radiation field. For relativistic electrons, typically, 7 >> 1 and k is 
nearly along the direction Z in Fig. 1., which is the direction into which a ray along if* 
would be reflected if the screen were replaced by a mirror. The vector k — k x + k y + k z 
k(9 x x + 9 y y + z) (see the x, y,z coordinate system in Fig. 1.), where 9 2 = d 2 x + Q 2 y and 
O x fiy « 1 . For \1> = 45° , we have k x m k(8 x — 7~ 2 /2) k9 x . Therefore, in general, the 
radiation pattern is shifted by an angle 7 _2 /2 << 1 . 
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The x and y integrations are straight forward. The variable y ranges from — oo to a\ and 
a 2 to +00, where the electron trajectory along 1? is taken to be a distance a\ from the lower 
edge and a distance a 2 from the upper edge of the slit so that a = a± + a 2 in Fig. 1. The x 
integration yields 2ir5{k' x — k x ) so that the k' x integral can be done trivially. Using k x fa k x , 
when 7 >> 1 , we obtain 

EM = ±l rm {L [h + , 2 ] + ( ^ + c>2) } (19, 
C tu ^-2^ I e-^ )a \ (20) 

-00 \ K x ' K y ' a ) K y K V 



where 



and 

/+00 7 1 

K lWTWT^y I ^~ ky)a2 - (21) 
■00 \ K x ' K y ' a ) K y K V 

Now we can write 

kl + k 2 y + a 2 = (k' y -tf)(k' y + tf), (22) 

where 

f = k 2 x + a 2 . (23) 

We see that ii and I 2 have poles at k' y = ±if ,and at k' y = k y , where ii = I%{o2 — * Qi)- 
These integrals can be done using the contours shown in Fig. 14. The results are 

'1 = - I^T^T-T' " 7 - „, .A , (24) 



/ (/ - ifc,,) (Ag + fc. 



2 + 



and 

Substituting these expressions for Ji and 7 2 into Eq. (19) above, one obtains Eqs.(2). Eq.(3), 
the expression for E y (k x , k y ) can be derived in a similar manner. Note that when a 12 — > 0, 
Eq. (19) gives the result for the transition radiation field as expected. 
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FigureCaptions 



Figure 1. A side view(lA), and a top view (IB), of backward diffraction radiation emitted 
by a charge passing through a slit which is oriented with its edge parallel to the plane 
of incidence (n, v plane). 

Figure 2. A side view(2A), and a top view (2B), of backward diffraction radiation emitted 
by a charge passing through a slit which is oriented with its edge perpendicular to the 
plane of incidence (n, v plane). 

Figure 3. Horizontal (3A) and vertical components (3B) of the intensity of DR for R = 
a/jX = 0.5; the intensity (Z) axis is scaled in proportion to the corresponding intensity 
component of TR. 

Figure 4. Horizontal (4A) and vertical components (4B) of the intensity of DR for R = 2.0; 
the intensity axis is scaled in proportion to the corresponding intensity component of 
TR. 

Figure 5. Effect of the divergence 8' on the horizontal component of the intensity. Shown 
are line scans in the variable X for two values of the variable Y = 8' = 0.2 and 
Y = 8 = 0, for a fixed value of the divergence e' — 0.1 and R = 2.0. 

Figure 6. Effect of divergence e' on the horizontal intensity. Shown are line scans in the 
variable X for a fixed value of the beam size 8 = 300/x and R = 2. 

Figure 7. Effect of beam size 8 on the horizontal intensity. Shown are line scans taken in 
the X, Z plane, i.e. Y — 0, for a fixed value of divergence e' — and R = 2.0. 

Figure 8. Effect of the beam divergence e' on the vertical intensity. Shown are line scans 
in the variable Y for two values of the variable X = e' = 0.2 and X = e' = 0, for a 
fixed value of the divergence 8' — 0.1, and R = 2. 

Figure 9. Effect of divergence 8' on the vertical intensity. Shown are line scans taken at 
X = for a fixed value of the beam size, 5 =200/i and R = 2.0. 

Figure 10. Effect of beam size 8 on the vertical intensity. Shown are line scans taken at 
X = for a fixed value of the divergence, 8' = and R = 2.0. 

Figure 11. Configuration of a DR interferometer composed of two slits oriented with the 
edges parallel to the plane of incidence. 

Figurel2. Vertical component of intensity from a DR interferometer for an electron beam 
with energy E = 500 MeV; R = 2, A = 3.2/im, 8 = 200/i, 8' = 0.05 and separation 
distance L = 6L V = 3 meters. 

Figurel3. Horizontal component of intensity from a DR interferometer for an electron 
beam with same parameters as in Fig. 12. 
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Figure 14. Contours used to evaluate the integrals I\ and I 2 denned by Eqs. (20) and 
(21), respectively. mo 
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